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^ ' Abstract 

c/3 ; 

An explicit expression is derived for the distribution function of end-to-end vectors and 
Cu ■ for the mean square end-to-end distance of a flexible chain with excluded- volume interactions. 

C I The Hamiltonian for a flexible chain with weak intra-chain interactions is determined by two 

small parameters: the ratio e of the energy of interaction between segments (within a sphere 
whose radius coincides with the cut-off length for the potential) to the thermal energy, and the 
Q ■ ratio S of the cut-off length to the radius of gyration for a Gaussian chain. Unlike conventional 

O I approaches grounded on the mean- field evaluation of the end-to-end distance, the Green function 

is found explicitly (in the first approximation with respect to e). It is demonstrated that (i) the 
distribution function depends on e in a regular way, while its dependence on d is singular, and 
(ii) the leading term in the expression for the mean square end-to-end distance linearly grows 
with e and remains independent of 6. 
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1 Introduction 

This paper is concerned with the distribution function of end-to-end vectors for a flexible chain 
with weak excluded-volume interactions. The effect of intra-chain interactions on the characteristic 
size of "real" flexible chains has attracted substantial attention from the beginning of the 50s of 
the past century, see early studies ^-0- A noticeable progress has been reached in this area 
during the past five decades, see |H1-I21], to mention a few works, and the results of investigation 
have been summarized in a number of monographs |25j-|27j. Nevertheless, the effect of segment 
interactions on the statistics of macromolecules has remained a subject of debate in the past five 
years. This conclusion is confirmed by a large number of publications, where this phenomenon 
is discussed, see J28j-j41|. and it may be ascribed to the importance of steric interactions for the 
description of such phenomena as (i) protein folding and denaturation 13 7 j , (ii) unzipping of DNA 
molecules |4()| . (iii) force-stretch relations for long chains at relatively small elongations j34j . and 
(iv) transport of flexible chains through narrow pores [231 • The account for intra-chain interactions 
(i) provides a consistent description of the dependence of intrinsic viscosity on molecular weight in 
dilute polymer solutions ^lElIlll) (h) leads to an adequate characterization of time-dependent 
shear stresses in polymer fluids [201 |^ ESI OOl UHl HEl IHHIi (hi) allows transport coefficients in 
polymer melts to be correctly estimated |33| . and (iv) predicts reinforcements of interfaces between 
incompatible solvents by random copolymers [SH]. Interactions between segments are assumed to 
be responsible for (i) swelling of polymer gels ^1 [211 , (ii) mechanical jSHl IHHj and optical ^21 
anisotropy of polymer solutions, (iii) structure of polyelectrolyte chains ^JElllSj, (iv) stability 
of globular proteins [2^1127], etc. 

There are two conventional ways to describe excluded-volume interactions in a fiexible chain. 
According to the first, a freely jointed chain is thought of as a set of rigid segments with given 
length linked in sequel, whereas the intra-chain interactions do not permit each pair of segments 
to occupy the same place. This model implies that the distribution function for end-to-end vectors 
coincides with that for a self-avoiding random walk. Although an explicit expression for the latter 
function is unknown, it is believed that the main term in the expansion of the distribution function 
for a long chain (with the number of segments A^ ^ 1) is described by either a stretched exponential 
function or a stretched exponential function multiplied by some power of the end-to-end distance 
j431 I44j . Results of numerical simulation confirm these hypotheses, but appropriate exponents are 
not precisely known |45j . 

According to the other method, a chain is treated as a sequence of beads bridged by (entropic) 
elastic springs [2S] . The Hamiltonian of the system equals the sum of the conventional Hamiltonian 
for a Gaussian chain and the excluded-volume functional, where the energy of repulsive interactions 
between segments is approximated by the Dirac delta-function. To assess the effect of steric inter- 
actions on the mean square end-to-end distance, it is presumed that the intensity of interactions is 
small compared to thermal energy. 

By applying various approximate techniques, the average end-to-end distance was evaluated in a 
number of studies (to the best of our knowledge, only a few attempts were undertaken to determine 
the distribution function of end-to-end vectors |H1 EI ) • It is widely accepted J2S1 |2H] that (i) the 
mean square end-to-end distance B of a flexible chain with segment interactions is determined by 
some small dimensionless parameter e (which quantity coincides with the product eS introduced 
later), (ii) the ratio (B/b)'^, where b is the mean square end-to-end distance for a Gaussian chain, 
may be expanded into a Taylor series in e, and (iii) the leading term (proportional to e) in this 
expansion is positive, which means that excluded-volume interactions are always repulsive. The 
common feature of the mathematical methods applied to develop these results was that all of them 
were grounded on the mean-field approach (the propagator for an appropriate path integral was 



evaluated on the classical path only). 

The objective of the present study is to develop an explicit expression for the distribution 
function of end-to-end vectors by evaluating the functional integral on (more or less) arbitrary paths 
(conformations of a chain). Our main result is that conclusions (i) to (iii) appear to be questionable. 
It is demonstrated that the Green function for a flexible chain with excluded-volume interactions 
is characterized by two small parameters, one of which, e, describes the intensity of intra-chain 
interactions per small volume associated with the segment length, while the other, 6, determines the 
ratio of the segment length to the average end-to-end distance for a Gaussian chain. The difference 
between these parameters is that the dependence of the distribution function on e is regular, 
whereas its dependence on 6 is singular. It is found that the leading terms in the expression for the 
mean square end-to-end distance is proportional to e, the scaling that substantially differs from the 
results of previous works. The growth of the average size of a chain with the intensity of segment 
interactions is ensured when the excluded-volume potential is included into the Hamiltonian with 
the sign different from the conventional one. 

The exposition is organized as follows. A Hamiltonian for a flexible chain with excluded- 
volume interactions is discussed in Section 2. This Hamiltonian is calculated on the classical path 
in Section 3. Mean- field approximations for the Green function and the end-to-end distance are 
given in Section 4. A perturbed Hamiltonian is determined in Section 5. Our analysis is based 
on the conventional approach, according to which the excluded-volume potential is expanded in a 
Taylor series in the vicinity of the classical path up to quadratic terms (a similar technique was 
recently employed in [IHI for the screened Coulomb potential). A leading term (with respect to 
small parameters) in the expression for the Green function is determined in Section 6. The mean 
square end-to-end distance B is evaluated in Section 7 for weak (e5^^ <^ 1) and in Section 8 
for moderately strong (arbitrary values of this ratio) interactions. Some concluding remarks are 
formulated in Section 9. Mathematical derivations are presented in Appendices A to F. 

2 Formulation of the problem 

A flexible chain is treated as a curve with "length" L in a three-dimensional space. An arbitrary 
configuration of the chain is determined by the function r(s), where r denotes radius vector and 
s G [0,L]. For definiteness, we assume the end s = to be fixed at the origin, r(0) = 0, and the 
end s = L to be free. An "internal structure" of the curve is characterized by a segment length bo 
and the number of segments A^ ^ 1, which are connected with length L by the formula L = b^N. 
In the absence of segment interactions, the Hamiltonian 

is ascribed to the chain. Here fee is Boltmann's constant, and T is the absolute temperature. The 
energy of intra-chain interactions is described by the functional 

V{r) = ^£ds £ V{r{s) - r{s'))ds', (2) 

where V{r) is a dimensionless potential energy of interactions, and vq is their intensity (which is 
assumed to be positive for repulsive and negative for attractive interactions). We confine ourselves 
to isotropic functions V^(r) = 14 (r) with r = |r|. The entire Hamiltonian of the chain reads 

H = Ho + V, (3) 



which iniphes that 



^(r) = k^T 
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V(r{s)-r{s'))ds' 



(4) 



Although Eq. (3) is employed in a number of studies, see |2S] and the references therein, the sign 
of the second term in this formula seems questionable, because V may be considered as a work of 
some fictitious external field that prevents segments from being located at the same positions. The 
latter implies that the contribution of V into H should be negative (as, for example, that of the 
work of an external force at extension of a chain). We postpone a discussion of the sign of V to 
Section 7, and suppose at this stage that vq may have an arbitrary sign. 
The Green function (propagator) for a chain G(Q) is given by 



G(Q) 



r(L)=Q 
r(0)=0 



exp 



H{r) 



ksT 



P(r(s)), 



where the path integral is calculated over all curves r(s) that obey the boundary conditions 

r(0) = 0, r(L) = Q. 



(5) 



(6) 



The aim of this study is to develop an explicit expression for the Green function of a chain with 
relatively weak segment interactions, when the last term in Eq. (4) is "small" compared to the 
first. The derivation is performed for an arbitrary function V{r), but our specific interest is in the 
excluded- volume potential V{r) = 5D(r), where 5Y){r) denotes the Dirac delta- function, and (5D(r) 
stands for its regularization after cut-off at small distances. To provide a rigorous definition, we 
introduce the Fourier transform of V{r), 



U{k) = fv{r)ex.p{ik-r)dr 



with the inverse 



V{r) 



(2vr) 



-— ^7 / f/(k) exp(— zk • r)dk, 
nr J 



(7) 



(8) 



where the dot stands for inner product. Equation (7) implies that for an isotropic function V{r), 
the function U(k) is isotropic as well, 

U{k) = U^{k), U^{k) = — / K(r) sm{kr)rdr. 

k Jo 

The Fourier transform of the Dirac delta- function 5D(r) reads ?7d*(^) = 1. To avoid divergence of 
the integral in Eq. (4), this function is cut off at an appropriate length-scale 6* = A;~^, 



UMk) 






(9) 



and the potential 6j){r) is defined as the inverse Fourier transform of UD^{k). To ensure that the 
potential ^/d* describes intra-chain interactions properly, we assume the length-scale of cut-off 6* 
to be substantially lower than the mean-square end-to-end distance for a Gaussian chain b = ^/h^L, 



6 = — <^l. 





(10) 



3 Weak excluded-volume interactions 



To characterize the "smahness" of the potential V, we, first, introduce a test function rg (the 
classical path) that minimizes functional (1) under conditions (6). Then the values of Hq and V are 
determined on the curve ro(s), and the dimensionless parameter e is chosen from the condition that 
y(ro) is small compared with HqItq) for any properly scaled end-to-end vector Q. We introduce a 
Cartesian coordinate frame {x,y,z}, whose z axis is directed along the vector Q, and present the 
minimization problem for the non-perturbed Hamiltonian Hq as follows: 
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dz . , 



ds, 



x{0) = y{0) = z{0) = 0, x{L) = Q, y{L) = z{L) = 0. 
The unique solution of this problem is given by 

ro(s) 
It follows from Eqs. (1) and (11) that 



Q-. 



Hoiro) = -kBTQ\ 



^-l 



Equation (8) implies that 

V{v{s) - r(s'; 



1 



(27r)3 



[/(k) exp — ?k • ( r(s) — r(s' 



dk. 



Substitution of this expression into Eq. (2) results in 



V{v) 



2L2(27r)3 
Combining Eqs. (11) and (13), we find that 

VQ 
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Simple algebra implies that (Appendix A) 
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U{'k)dk I ds f exp(^-tk-Q^-^]ds'. 



^(ro) 



2(7r6Q)2 Jo x^ Jo 



1 dx '■°° 
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(12) 



(13) 



(14) 



(15) 



Formula (15) determines the energy of intra-chain interactions for an arbitrary function U^{k). In 
particular, for the excluded-volume potential (9), this equality reads 



^(ro) 
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Calculation of the integral (Appendix A) results in 



/ 4 / \'i~-cos{kbQx))dk. 



(16) 



y(ro) = ^Aik^bQ), 



(17) 



where 

A{x) = \ r{x^ - z^)'^-^^dz. (18) 

x'^ Jo z^ 

The function A(x) is plotted in Figure 1, which shows that A monotonically decreases with x. The 
hmits of A{x) as x ^ and x ^ oo are given by (Appendix A) 

hm A(x) = -, hm A(x) = 0. (19) 

It follows from Eqs. (12) and (17) that the smallness of the functional V compared with the 
non-perturbed Hamiltonian Hq is tantamount to that of the pre-factor 

vok^ vo 



compared with thermal energy ksT. Here v = ovrfo^ is the volume in which interactions between 
segments are taken into account. Introducing the parameter e by 

- "° (20) 



UbTv 



we conclude that the condition |e| ^ 1 means that the energy of excluded- volume interactions at 
the length-scale of cut-off is small compared with thermal energy. Substitution of expressions (12), 
(17) and (20) into Eq. (4) results in 



i^(ro) = teT[|Q^ + ^A(|; 



(21) 



According to Eq. (21), in the mean-field approximation, the functional H is determined by two 
dimensionless quantities, e and 5. The former characterizes the smallness of excluded- volume inter- 
actions, whereas the latter describes the length-scale where these interactions are important. The 
parameter e is located in the numerator, which means that it is responsible for regular perturbations 
of Hq, while 5 stands in the denominator, which implies that it describes singular perturbations. 
This distinguishes the present approach from previous studies, where one small parameter was 
introduced, and its effect was assumed to be regular in the sense that the Green function was 
expanded into a Taylor series with respect to this quantity. 

Before proceeding with the calculation of the Hamiltonian H on an arbitrary curve r(s), it is 
instructive to evaluate the Green function G(Q) and its second moment on the classical path ro(s). 

4 Mean-field approximation of the Green function 

Our aim is to find an approximation of the mean square end-to-end distance B when the set of 
admissible curves {r(s)} contains the classical path ro(s) only. Substitution of expression (21) into 
Eq. (5) results in 



G(Q) = Cexp[-(^Q2^|-^(|: 



(22) 



where the pre-factor C characterizes the measure P(r(s)) of the function ro(s). To determine C, 
we treat the Green function as the distribution function of end-to-end vectors Q that obeys the 
normalization condition 

j G{Cl)dCl = l. 



In the spherical coordinate frame {Q,(j),9} whose z vector is directed along the vector Q, this 
equality reads 

("OO 

47r / G{Q)Q'^dQ = 1. (23) 

Jo 

Substituting expression (22) into Eq. (23) and setting Q = bQ, we find that 

'3Q2 £ Q XI ~2 ~ 



Anb'^C / exp 
Jo 



Expanding the function under the integral into the Taylor series (this is possible because the 
function A(x) is uniformly bounded) and disregarding terms beyond the first order of smallness 
with respect to e, we obtain 



47r63cy'"exp(-^ 



37r ^6' 



Q^dQ = 1. 



It follows from this equality that with the required level of accuracy, 



where 



C. = l Aif)e.p( ^ 



The leading term in the expression for Ci reads (Appendix B) 

7r(5 



Ci 
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(24) 
(25) 

(26) 



The mean square end-to-end distance of a chain is given by 

B"^ = fq^GiCDdQ. 
In the spherical coordinates {Q,(f),9}, this equality is presented in the form 



B^ = 4Tr G{Q)Q''dQ. 
Jo 



(27) 



Substituting expression (22) into Eq. (27), introducing the new variable Q = Q/b, and neglecting 
terms beyond the first order of smallness with respect to e, we find that 



B^ = ATTb^C 



' 3Q2 
exp' 



/o ^ 2 

This equality implies that with the required level of accuracy, 



\-^A(^)\Q^dQ. 



(|)^.c(^)i(,_J4B, 



where 



B, 



^(|)exp(-^)Q4dQ. 



(28) 
(29) 



The leading term in the expression for Bi reads (Appendix B) 



B. = f. (30) 

It follows from Eqs. (24), (26), (28) and (30) that 

(f)'=l + .y4(C,-B,) = l + ^/f, (31) 

V / V TT'^ 9 V 27r 

Equation (31) implies that at vq > 0, the mean square end-to-end distance of a flexible chain 
with excluded-volume interactions exceeds that of a Gaussian chain. At first sight, this conclusion 
confirms that the sign of the contribution of V into the Hamiltonian H is chosen correctly in Eq. 
(3). Setting 6* = bo (the cut-off of the potential C/* occurs at the segment length) and using Eqs. 
(10) and (20), we find that 

(^Y = l + ^° , 

^b ^ ^/OTrkBTutuhc 

where ftubc = 47r6ofe denotes volume of the characteristic tube around a chain (a circular cylinder 
whose radius coincides with the segment length and whose length equals the mean square end-to- 
end distance). Our aim now is to demonstrate that Eq. (31) does not capture the main contribution 
into the asymptotic expression for the average end-to-end distance. To prove this fact, we derive 
an explicit expression for the Hamiltonian H that accounts for second order terms with respect to 
admissible fluctuations from the classical path ro(s). 

5 Perturbations of the Hamiltonian 

In accord with Eq. (6), the function r(s) is given by 

r(s) = Q| + R(s), (32) 

where the function R(s) satisfies the boundary conditions 

R(0) = 0, R(L) = 0. (33) 

To simplify the analysis, we neglect longitudinal fluctuations and expand the function R(s) that 
describes transverse fluctuations into the Fourier series 

oo oo 

r^(*) = 2^ Xm sm ^— ei + 2^Ym sm ^—62, (34) 

m=l m=l 

where Xm, Ym are arbitrary coefficients, and ey^ [k = 1,2,3) are unit vectors of a Cartesian 
coordinate frame, whose 63 vector is directed along the vector Q. Any function R(s) obeying Eq. 
(34) satisfies also boundary conditions (33). We substitute expressions (32) and (34) into Eqs. (1) 
and (2) and, after some algebra, find that (Appendix C) 

i?o(r) = i/o(ro) + ^^E"^'(^™ + ^™)' 

m=l 

00 

V{v) = y(ro)--^ Y. PmniX^Xn+Y^Yn), (35) 

^ ' m,n=l 



where 
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X F{kQx + 7rn) — F{kQx — irn 
+F{kQx — Trm) F{kQx + vrn) + F{kQx — irn 
and the function F{z) reads 



F{kQx + TTTTl) 



F(^) 



1 z 
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(36) 



(37) 



It follows from Eqs. (4) and (35) that with the required level of accuracy, the perturbed Hamiltonian 
H is given by 



H{r) = F(ro) + k^T 



3lT 



2 oo 



J2^'(Xi + Y, 
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m=l 



^ ' m,n=l 



(38) 



Our purpose now is to substitute this expression into Eq. (5) and to calculate the path integral. 



6 The Green function 

Combining Eqs. (5) and (38), we find that 
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P(R(s)). 



(39) 



Bearing in mind Eq. (21) and using the matrix presentation of the functional integral, we obtain 



G(Q) = C7exp 



3Q 



,Q 



2 +3^^^T^)J'^'^^)' ^^^^= ^""^^ 
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dX, 



where C is a constant associated with transition from the measure P(R(s)) to the measure dX 
(this quantity will be found from the normalization condition for the Green function). Here X is 
the vector with components X^, A is the diagonal matrix with components Amm = 37r2r7T,2/(262), 
and B is the matrix with components Bmn = —'^Pmn/i^Tr'^). Calculation of the Gaussian integral 
implies that 

^3(52 e ,,<5^M detA 



G(Q) = C7exp 
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In the first approximation with respect to e, 
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It follows from Eqs. (40) and (41) that 
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Substitution of expression (36) into Eq. (43) implies that (Appendix D) 
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The infinite sums are calculated in Appendix D, where it is shown that 



'S'o(^) = ^"(^)> 



, , smz 
a{z) = 1 + 2 6- 



(45) 
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It follows from Eq. (46) that Sq{z) is an even function. Using this property, we present Eq. (44) 
in the form 



vr 



2 rl 



l-x^ 



Q2 7o x2 



dx / U^{k)a{kQx)k dk. 



(47) 



Formula (47) determines the denominator in Eq. (42) for an arbitrary isotropic potential V^(r). 
For the function f^D*(fc) given by Eq. (9), this equality reads (Appendix D) 



«-^-o(f). 



where 



Mx)=-j_{x--z-r^'^ 



1 

x^ Jo 



dz. 
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The function ^o(^) is plotted in Figure 1. This function is even, it is negative for any x £ (—00, 00), 
and it monotonically increases with \x\. The limits of the function Ao{x) are given by (Appendix 

D) 

2 
lim Ao{x) = - — , lim Ao{x) = 0. (50) 



x->0 
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It follows from Eqs. (42) and (48) that 

G(Q)=Cexp[-(^ + ^A( 



3Q!^^ Q. 
2 Svr ^6 



^ + 36^^:^" 



(51) 
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Equation (51) provides an analytical expression (up to the normalization constant C) for the Green 
function of a flexible chain with excluded- volume interactions. The only assumption employed in 
the derivation of this formula is the smallness of e compared with unity (no limitations on 5 were 
imposed). Given a coefficient (5, Eq. (51) is valid provided that the expression in the last square 
brackets does not vanish. 

Our purpose now is to determine the normalization constant C and the mean square end-to-end 
distance B. We begin with the case of "weak" interactions when the dimensionless parameter e5~'^ 
(and, as a consequence, e) is small compared with unity. 

7 Weak segment interactions 

To determine the normalization constant C, we substitute expression (51) into Eq. (23) and obtain 



A^Tb^^C I exp 
/o 



/3Q2 



Q- 



V 2 +3^^(T 



^ + 36^<^" 



Q^dQ = 1. 



Expanding the function under the integral into the Taylor series and neglecting terms beyond the 
first order of smallness with respect to e and £(5^^, we find that 



Airb'^C I exp 
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3Q^ 
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It follows from this equality that with the required level of accuracy. 
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where Ci is given by Eq. (25), and 

C2 = 



Ao(|)exp(-^)Q2dg. 



The leading term in the expression for C2 reads (see Appendix E for detail) 

Co = — \ — S . 
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Substituting expressions (26) and (54) into Eq. (52) and neglecting terms beyond the first order of 
smallness, we arrive at the formula 
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To calculate the mean square end-to-end distance B, we substitute expression (51) into Eq. (27), 
set Q = Q/b, disregard terms beyond the first order of smallness with respect to e and e6~'^, and 
find that 



B^ = A-kWC I exp 
/o 



It follows from this equality that 
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where Bi is given by Eq. (29), and 

B2 = f ^o(|) exp{-^)Q^dQ. (57) 

The leading term in the expression for B2 reads (Appendix E) 

4 [2^ 



Substituting expressions (30), (55) and (58) into Eq. (56) and neglecting terms beyond the first 
order of smallness, we find that 

With reference to the conventional standpoint [the positive contribution of V into the Hamiltonian 
H in Eq. (3)] Eq. (59) contradicts our intuition: repulsive interactions between segments result in 
a decrease in the end-to-end distance. Formula (59) contradicts also the mean-field Eq. (31): the 
ratio [B/b]'^ is proportional to e instead of the classical scaling e6. 

The latter discrepancy may be explained if we recall that excluded-volume interactions do not 
permit different segments of a flexible chain to occupy the same positions, which means that their 
effect is substantial only for "curved" configurations of a chain, whereas the mean-field approach is 
confined to "straight" configurations. This implies that the mean-field technique is inapplicable to 
the analysis of a fiexible chain with long-range interactions between segments, in agreement with 
the conclusion derived about 50 years ago |2] (based on different arguments). 

Equation (59) leads to a physically plausible result (an increase in the average end-to-end dis- 
tance with intensity of excluded-volume interactions) when e is negative. The latter is tantamount 
to the negativity of contribution of V into the Hamiltonian H in Eq. (3). To demonstrate the 
correctness of this assertion, it is instructive to compare the values of H on two paths: (PI) a 
straight line (11) that connects the end-points, and (P2) a sinusoidal path (32) and (34) with the 
only non-zero term corresponding to a fixed m > 1. According to the physical meaning of excluded- 
volume interactions, the energy H2 on path (P2) should exceed the energy Hi on the straight path 
(PI) as Q ^ 0, i.e., when the chain is superposed on itself several times. On the other hand, it 
follows from Eq. (35) that 



H2 = Hi + 



— 7I2 TtV / UD*{k)k'^dk / (1 - x^)pmm{kQx)dx 



{Xl + Yl), (60) 



where Pmm is given by Eq. (C-29). The integral in Eq. (60) can be calculated explicitly, but this 
is not necessary for our purpose, because Eq. (C-29) implies that in the limit of at small Q, the 
function Pmm reads 

Pmm\kQx) ~ L . 

Combining this estimate with Eq. (60), we see that Eq. (3) with the positive contribution of V 
results in H2 < Hi for sufficiently large vq, which contradicts the definition of excluded-volume 
interactions. 

Based on this analysis, we conclude that the potential of intra-chain interactions should be 
included into Eq. (3) with the negative sign, and, as a consequence, the distribution function of 
end-to-end vectors (51) should read 



^Q 2 ,2^^Q^]l. . rQ 



2 3"^ '6 
12 



G(Q) = Cexp -^ + ^/x52^(^) l + f^Ai{^) \ (61) 



S 



-1 



where we introduce the notation 



M 



27r(52' 



Ai{x) 



--,M-) 



(62) 



and preserve the positiveness of e ^ 1. We intend now to derive exphcit expressions for the nor- 
maUzation constant C and the mean square end-to-end distance B for an arbitrary (non necessary 
small) values of the dimensionless parameter n. Our aim is to demonstrate that Eqs. (55) and (59) 
provide the leading terms in the expressions for C and B (after an appropriate corrections of signs) 
when the condition e5~'^ <C 1 is violated. 



8 Moderately strong interactions 

It follows from Eqs. (23) and (61) that 



Bearing in mind that 



Airb'^C I exp 
/o 



1 + iM'^d : 



3x2^1+1^^^(1) 

— ^—x ax = 1. 



2 J 1 + Mi(f) 



1 + Mi(f] 



l-fi 



A,if)-l5'A{f] 
1 + Mi(f) 



/ 3 xl 



we find that 

^=l2^J 
Evaluating the integrals 

Ml 



1-6^ 



Ai 



:f ) - ¥'Mj] 
i+Mi(f) 



3x' 



X dx 



Ai 



l + /uAi(f) 



3x^ 



X dx, 



Ao 



Ml) 



1 + Mi(f) 



3x' 



x'^dx, 



and neglecting small terms, we arrive at the formula (see Appendix F for detail) 

"2 / 3 \|e /""^ Al{x)x'^ 



^"(2^ 



^'l + i 

vr 



^_/£\l£ 



_9 V27r/ 5 Jo l + fiAi{x) 
Substituting Eq. (61) into Eq. (27), we find that 



dx 



B\2 
b 



(!) 



2vr6M^ 






1 + M^i(f) 



3. -Vg7 ^_^^^^4^^ 



The integrals 
El 



°° ^i(f) 
1 + Mi(f) 



3x^ 



x dx, 



^(f) 



i + Mi(f) 



3x^ 



X dx 



(63) 



(64) 



(65) 



(66) 



(67) 



are calculated in Appendix E, where it is shown that the leading term in the expression for the 
mean-square end-to-end distance B is given by 



(fr^(^)H-i 



81 FYe'^ '■°° 



AK 



XjX 



87r2 V 27r 5 Jo I + fJ.Ai{x 



-dx 



(68) 
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Equations (65) and (68) imply that the leading term in the expression for the mean square end-to- 
end distance reads 



B\2 e 



Q 



/ TT 



(69) 



.27 V27r/ 5 Jo l + ^iAi{x) / 

Formulas (65) and (69) differ from Eqs. (55) and (59), respectively, by the integral terms in the 
square brackets only (after correction of the sign of e). Although these terms contain the pre- 
factor e/5 which may accept arbitrary values, simple algebra demonstrates (Appendix F) that 
their contributions are negligible. 

9 Concluding remarks 

The formula 



2e- 



3Q2 e ,,Q,ir, e ^ ,Q 



has been derived for the distribution function of end-to-end vectors for a flexible chain with weak 

excluded-volume interactions. Here b is the mean square end-to-end distance for a Gaussian chain, 

e is a small parameter that describes the intensity of segment interactions, 5 is the ratio of the 

average end-to-end distance for a Gaussian chain to its segment length, and the functions A{x) and 

Ai{x) are given by Eqs. (62), (A-8) and (D-24). The leading term in the expression for the mean 

square end-to-end distance reads 

'B\2 4e , , 

Equation (70) differs from similar relations developed in previous studies, where the ratio on the 
left-hand side was found to be proportional to eb. A reason for this difference is that our result is 
grounded on the calculation of an appropriate path integral, whereas conventional conclusions are 
obtained by using mean-field approximations. It appears that the latter approach is inapplicable to 
problems where interactions between segments located far away (along a chain) from each other are 
substantial, because the mean-field technique is confined to "straight" paths only, while inter-chain 
interactions reveal themselves mainly on "curved" configurations of a chain. 

It has been found in the calculation of the path integral that the sign of the excluded-volume 
potential should be corrected in the conventional formula (3) for the Hamiltonian. It has also been 
shown that the mean square end-to-end distance cannot be expanded into a Taylor series with 
respect to small parameters (as it is traditionally presumed), because even sub-leading terms of the 
highest order in the formula for B^ include (5 In 5, see Eq. (E-11), and \/e, see Eq. (F-25). 

This work focuses on the analysis of the distribution function for a flexible chain with excluded- 
volume potential (9). However, the results can be easily extended to an arbitrary potential of 
intra-chain interactions by using Eqs. (15) and (47). In particular, these relations allow an explicit 
formula to be derived for the distribution function of end-to-end vectors for flexible polyelectrolyte 
chains. The latter will be the subject of a subsequent publication. 
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Appendix A 

To transform integral (14), we choose a spherical coordinate frame {k, (f>, 9}, whose z axis is directed 
along the vector Q, and obtain 

Vn f°^ r"^^ f^ f^ f^ f S — s' \ 

Vivo) = / U^:(k)k dk dcp sinOdO / ds ex.p(—tkQcos6 — - — ]ds'. 

2L^[2Tr)-^ Jo Jo Jo Jo Jo ^ ^ ^ 

Calculating the integral over (f) and introducing the new variable x = cos 9, we find that 

F(ro) = ^^J^ y°° U,{k)k^dkJ J{kQx)dx, (A-1) 

where 

J{z) = ds expf — iz — - — jds' = / expf — j~)ds / expf — — jds'. (A-2) 

For an arbitrary non- negative z, we have 

Jo exp(-^)cis = -(l-expHz)), J^ exp(^)ds' = -(exp(.z) - l), (A-3) 

It follows from Eqs. (A-2) and (A-3) that 

J(z) = ^-(l- cos z). (A-4) 

z^ 



Substitution of Eq. (A-4) into Eq. (A-1) implies (after changing the order of integration) that 

^^"■"^ = 2^ to % r ^* ^^^ (^ " ^°'(^^^)) ^^- ^^-^^ 

Combining Eqs. (12) and (A-5), we arrive at Eq. (15). 

Introducing the new variable z = kbQx, we find from Eq. (16) that 



vo /■! dx r^'bQx^ 
27r2(6g)3 Jo ^ Jo 



^i^o) = ^ ^,,^,^ L ~ I (1 ~ cos z)dz. 



Setting y = k^bQx, we obtain 



2TT^bQ Jo y^ Jo 



V(ro) = :P^I ^^ I (1- COS z)dz. 



We now change the order of integration 



/ (1 — coszjdz / -^ 

2tt%QJo Jz y^ 



y(ro) = ^P^ I (1 - cos z)dz I ^, (A-6) 



and calculate the internal integral 

rk,hQ ^y {k^hQf - z2 



y3 2z2(fe,6Q)2 • 
Substitution of this expression into Eq. (A-6) implies Eqs. (17) and (18). 
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To determine the limits of the function A(x), we transform Eq. (18) as follows: 

1 f"^ 1 — COS z 1 f^ 

A(x) = - K dz ^/ (1- cos z)dz. (A-7) 

X Jo Z^ X'^ Jo 

The first integral is simplified by integration by parts, 



X 



1 — cos z , X — sin X 



-dz = 7^ h 2 / 5 dz 



z"^ x'^ Jo z^ 



Calculation of the second integral in Eq. (A-7) implies that 

rx 

/ (1 — cos2;)dz = X — sinx. 
Jo 

Substitution of these expressions into Eq. (A-7) yields 

^/ X 2 /"^ . , , , , z - sinz , ^ 

A[x) = — a(z)az, a[z) = ^ . (A- 

X Jo z-^ 

To find limits of the function A(x) as a; ^ and x -^ oo, we apply L'Hospital's rule 



2 1 — cos X 

- lim :; 

3 ^ x^ 



lim A{x) = - lim . (A-9) 



Equation (A-9) implies that 

1 .. sinx 1 



lim A{x) = — lim 



x^O 3 a;-^0 X 3' 

which coincides with the first equality in Eq. (19). The other equality follows from Eq. (A-9). 

Appendix B 

Equation (24) follows from Eq. (23) and the formula 



exp[-^)Q'dQ = -^ — . (B-1) 

To determine the constant Ci, we substitute Eq. (A-8) into Eq. (25), introduce the new variable 
X = Q/S, change the order of integration, and find that 



Ci = 26^ / a{z)dz / exp(--(5\^)xdx. 



a^z^dz I exp(- 
Calculating the internal integral. 



exp[--6^x')xdx = —exp[--6'z^), (B-2) 

and setting x = 6z, we arrive at 

2 r°^ / 3^2 2\ / A , 2 Z""^ ^x, / 3x2x 

3 / ^^PV~2 /'^^^'^ 3 / '^^^-''^^PV — 2~/ 

It follows from Eq. (A-8) that the function a(x) is even, which implies that 
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As the limit of the function a{x) as x ^ is finite and this function decreases being proportional 
to x^^ as X ^ oo, the Fourier transform 



a(s) 



a{x) exp(zsa;)(ix 



(B-4) 



exists, and 



— / a(s)exp(-— jds. (B-5) 

Substituting expression (B-5) into Eq. (B-3) and changing the order of integration, we obtain 



<> 



1 f°^ 
Ci = — / a{s)ds 



exp 



/3x^ isx 
■(- + - 



dx. 



To calculate the internal integral, we set y = x\f?> and find that 



exp 



/3x^ isx 



V 2 



b 



dx 



2tt / sS 



(B-6) 



(B-7) 



Substituting expression (B-7) into Eq. (B-6) and introducing the new variable x = 3/((5v3), we 
obtain 



Ci 



'i{6xV3) exp ( — —]dx. 



(B-8) 



3\/2tt j-oo 

If the function d{s) have had finite derivatives at s = 0, the integral were evaluated up to an 
arbitrary level of accuracy with respect to 5 by the stationary phase method. As this is not the 
case, we present Eq. (B-8) in the form 



3^/2^ 
calculate the first integral, 

and evaluate the residual 



— — ]dx+ (d((5xv3) — a(0)j expf — -— jdx 



Ci = -a(0)+i?i, 



(B-9) 



^1 



3V2TT 



a((5s v3) — a(0) ) exp ( ]ds. 



It follows from Eq. (B-4) that 

d{5sy/3) - d(0) = / a(x] 
As the function a(x) is even, this equality reads 



exp(?(5sxv3) — 1 



dx. 



d((5s\/3) — d(0) = / a(x) cos(5sxv3) — 1 



dx = —4 / a x)sm dx. 

Jo ^ ^ 2 



Substitution of Eq. (A-8) into this formula results in 



a((5s\/3) - a(0) = -36^s 



2„2 



sm 



SsxyS \ 2 



Ssxy/S 
2 



smx\ , 
1 ]dx. 



(B-10) 
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It follows from this equality that 



\aiSsV3) - a(0)| < 36 



h' 



sm 



Ssx\/3 \ 2 



5sx\/3 
2 



sinx 



dx 



j-C 



sm 



<5sa:\/3 \ 2 



5sx\/li 
2 



dx = 4(5|s|\/3 / ( -] dy = 2^5\s\V2,. (B-11) 



where we set y = \5sx\/?> and use the identity 



/■°^ /siny\2 vr 

Jo ^ y ^ 2 



(B-12) 



Combining Eqs. (B-10) and (B-11), we conclude that i?i is of order of 6^. Neglecting terms beyond 
the first order of smallness with respect to 5, we find from Eqs. (B-4) and (B-9) that 



2(5 f°° 
Ci = — / a{x)dx. 
3 Jo 

Integration by parts with the help of Eqs. (A-8) and (B-12) results in 



(B-13) 



a{x)dx = — . 
4 



Equations (B-13) and (B-14) result in Eq. (26). 

Equation (28) follows from Eq. (27) and the identity 



(B-14) 



/:-(-f)«'«4;!- 



(B-15) 



To determine the coefficient Si, we substitute expression (A-8) into Eq. (29), set x = Q/6, change 
the order of integration, and find that 



/•OO /'OO 

Bi = 26^ / a{z)dz / exp 

JO Jz 



3(5^x 



2^2, 



X dx. 



Bearing in mind that 



exp 



3,52x2 



x dx 



1/2 2 



3<52^2 



(B-16) 



we obtain 



4 r roo 
Bi = - 
9 



"(f) "^^-"1-)^^ + 2/0 «(|)e^p(-^jy'^y 



'y^ „-._,/ Va^_. , 3 r ^,y^_^( Sy\ 2 



(B-17) 



The first integral is given by Eq. (B-3). To determine the other integral, we use Eq. (B-5) for the 
function a and change the order of integration. 



.y 





y)e^p{-^)y'dy = - f a{s)ds 



2 r " AttJ-oo 
The internal integral in Eq. (B-18) is calculated explicitly 

'3y2 isy 



exp 



3y2 isy 



y'^dy. 



exp 



2 + 6 



2 , 1 /27r / s \/-, s 



(B-18) 



(B-19) 
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Substitution of Eqs. (B-3), (B-18) and (B-19) into Eq. (B-17) results in 



3 



6 f^^ / z^ \ 



(B-20) 



where z = s/{5\/i). Evaluating the integral by using the same approach that was applied to derive 
Eq. (B-13) and utilizing the identity 



exp 



(-'){l-z^)dz = Q, 



(B-21) 



we find that the order of smallness of the second term on the right-hand side of Eq. (B-20) is higher 
than that of the first term. Equation (30) follows from Eqs. (B-13), (B-14) and (B-20). 



Appendix C 

Substitution of expressions (32) to (34) into Eq. (1) results in 



Hniv) 



3kBT 



2bo [ L 



Q' 



+ E 



m,n=l 



-irm -nn 



■Kins irns 



{XrnXn + YmYn) / COS —— COS 

JO -t^ 



L 



ds 



Employing the orthogonality of the trigonometric functions and Eq. (12), we arrive at 



i/nfr) 



3/cbT r^2 ^ 



262 



2 oo 



Q' + ^E"^'(^™ + ^r 



m=l 



(C-1) 



The first equality in Eq. (35) follows from Eq. (C-1). 

Substituting Eq. (32) into the exponent in Eq. (13), expanding the obtained expression into 
the Taylor series, and omitting terms beyond the second order of smallness with respect to |R|, we 
find that 



exp 



ik Q 



s — s 



L 



-ik. • (r(s) — r(s') 

-i[k.(R(.)-R(.') 
Bearing in mind that Q = Qe^, we conclude from Eq. (34) that 



1 - ?k • (r(s) - R(s' 

2 



(C-2) 



k • Q = /csQ, k • R(s) = Y^ (kiX^ + k2Y^) sin 



m=l 



Trms 

IT' 



Substitution of these expressions into Eq. (C-2) results in 



exp 



ih ■ (r(s) — r{s'] 



exp 



i^-iks 



s — s 
) 

L 



oo 



1 - « ^ {kiXm + k2Yrr, 



X sm ■ 



L 



sm- 



L 



Y^ {kiXm + k2Ym) f sin — sin 



m,n=l 



m=l 

irms 



irms 



L 



L 



x(fciA„ + fc2i'^n)( sm^ sm 



L 



L 



(C-3) 
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Combining Eq. (C-3) with Eqs. (13) and (14), we find that 

V{r) = F(ro) + ^^2^ [^1 (R-) + ^2(R-)] , (C-4) 

where 

Vi= J U{k)Pi{k)dk, V2 = J U{k)P2{k)dk. (C-5) 

The functions Pi and P2 in Eq. (C-5) are given by 



00 



Pi(k) = --^{kiXm + k2Ym)pm{k:iQ), 



2 



m=l 

00 



P2{k) = -i ^ {klX^ + k2Yra){kiXn + k2Yn)pmn{kzQ), (C-6) 

m,n=l 

where the coefficients pm and pmn read 



PmiksQ) = ds expi^-iksQ — - — j (^sin — sin— —jiis', 

/, ^N f^ 7 f^ f , ^■5-s'\/ • '^"^■5 . vrms'x/ . vrns . vrns'x , , , „ _, 
Pm.n{kzQ) = j as exp\^-ik3Q — - — j (^sm — sm——j{sm— sm — — jds . (C-7) 

In spherical coordinates {k, (f), 6}, the Cartesian components of the vector k read 

ki = k sin 6 cos (p, k2 = k sin 9 sin (/>, k^ = k cos 0. 
Substitution of these relations into Eqs. (C-5) and (C-6) implies that 

Vi = -^Y1 I U^{k)k^dk I ^(x„cos(/. + y„sin0)# r pm{kQ cos 6) siii^ OdO = 0. (C-8) 

m=l 

By analogy with Eq. (C-8), we write 

t^2 = -- V / U,{k)kUk X^XnCOS^(t) 

-|-(Xm,yn -|- XnYm) sin (/) COS (j) + YmYn S\V? (f) d(f) I Pmn{kQ COS 9) sin^ 0d0. 

J JO 

Calculating the integral over (j) and introducing the new variable x = cos 9, we arrive at 

V2 = -- E {XmXn + Y„,Yn) U,{k)kUk Pmn{kQx){l-X^)dx. 

m,n=l J—L 

Substitution of expressions (C-8) and (C-9) into Eq. (C-4) implies that 



(C-9) 



CXj n r>C) p 1 

y(r) = y(ro)-:;:^^ 5] (X„X„ + y„y„) / U4k)kUk I pmn{kQx){l - x^)dx. (C-10) 

n,n=l ~ 



327r2L2 



Our aim now is to determine the coefficients Pmn- It follows from Eq. (C-7) that 

nj - „(!) _ ^(2) _ J2) , (3) /n 1 n 
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where 



p(i) 

fmn 



fmn 



/ exp(— i/csQ— ) sin — - — sin— — ds / explt/csQ— Ids', 
JO ^ -ty' L, L Jo V L/ 

/ expi^-tksQ- 



sm—— as / explzfcsy— sm — - — as, 
L Jq V L/ L 



p(3) 
fmn 



We begin with the quantity 



/ ex.p[—tksQ—]ds expfzfesQ— ) sin — - — sin — - — ds' . 
Jo ^ Ly Jo ^ Ly L L 

;ity 

[^ / , ^s \ Tims Tins l/rifn r,i2\\ 

Bmn = J expi^-tksQjj sm -j- sm -j^ds = -(^U - ^mn)> 



(C-12) 



(C-13) 



where 



mn 



B^H 



Using the formula cos ^ = 2 
we obtain 



expf— z/csQ— ) cos ds, 

^ L^ L 

= / exp -^fc3(5- cos— ^ — - — —ds. (C-14) 

Jo ^ Ly L 

expf ^ j + expf-^ j and calculating the integral in Eq. (C-14), 



bW 



2 Jo 
L 

~2i 



exp(-i{ksQ - 7r(m - n)) — ) + expl-i^ksQ + 7r(m - n)) — 



ds 



1 



-(ex.p(—i{ksQ — -iT{m — n)) 



+ 



ksQ + 7T{m — n) 
Taking into account that 



k^Q — 7r{m — n 

exp( —i{k^Q + 7r(m — n)) ) — 1 



1 



exp(— zx) — 1 = — 2zsin — exp(— — ), 



(C-15) 



we find that 



b'>^^ = Lexp[--k3Qj exp(^-i{m-n)jF{k3Q -TT{m-n)) 

+ exp( — —i{m — n)\F{kzQ + 7r(m — n)) 
where the function F{z) is determined by Eq. (37). Similar transformations result in 



5(2) 



^exp(--A:3Q 



Qxp{—i{m + n)]F{k3Q — 7r{m + n)) 



vr 



+ exp( i{m + n)]F{k3Q + ^{ra + n)) 



Substitution of these expressions into Eq. (C-13) implies that 



B„ 



-^exp(--A;3(5 



exp 



( —i{m — n) \ F^k-iQ — ii{m — n)) 



-|-exp( i{m — n)]F{k3Q + 7r{m — n)) 

exp(— z(m + n)jF{k3Q — T:{m + n)) 



TT 



+ expl-—i{m + n)\F{k3Q + 7r(m + n)) 



(C-16) 
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It follows from Eqs. (37), (A-3) and (C-15) that 

J exp(ik3Qj)ds' = 2Lexp{^k3Q)F{k3Q). 
Substitution of Eqs. (C-16) and (C-17) into Eq. (C-12) results in 



(C-17) 



„(i) 



L^FiksQ) 



expf — z(m — n)]F{k3Q — Tr{m — n)) 



+ exp( i{m — n)]F{ksQ + 7r(m — n)) 

expf — z(m + n)\F{k^Q — -ir{m + n)) 
+ expf-— z(m + n)jF{k3Q + 7r(m + n)) 
Replacing the exponents by trigonometric functions, we find that 



„(i) 



L'^F{k3Q){\ COS 



Tr{m — n) - 



FiksQ + 7r(m - n)) + Fik^Q - 7r(m - n)) 



cos ■ 



+?< — sin 



TT{m + n) - 



F{k^Q + 7r(m + n)) + F{k^Q - 7r(m + n)) 



7r(m — n) 



Fik^Q + 7r(m - n)) - ^(A^sQ - 7r(m - n)) 



+ sm ■ 



Tr(m + n) 



F{k^Q + Trim + n)) - F{k^Q - 7r(m + n)) 



We proceed with transformation of the quantity 



Cmn = / eyivy-kzQj) 



L/^^" — 



Tims' Tins' 
sin ■ 



L 



ds . 



Comparison of expressions (C-13) and (C-19) implies that Cmniks) coincides with B„ 
cording to Eq. (C-16), this means that 



(C-18) 



(C-19) 



-fe). Ac- 



On 



L 



■exp 



[n^sQjl exp(^-z(m- n)jF{-k3Q - 7r(m - n)) 



+ expf t(m — n)jF{—k^Q + Ti{ra — n)) 



vr 



expf — z(m -1- n)\F{—k^Q — 7r(m -|- n)) 



-fexp( i{m + n)]F{—ksQ + TT{m + n)) 

Bearing in mind that the function F{z) is even, we find that 



On 



-exp^-fcsQ 



exp(—i{m — n)\F[k^Q + 7r(m — n)) 

-|-exp( i{m — n)jF{ksQ — T:[m — n)) 

exp(-z(m + n)\F{k-iQ + 7r(m + n)) 



-|-exp(-— «(m + n)\F{k^Q - Tr{m + n)) 



(C-20) 
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It follows from Eq. (C-17) that 

rL 



J exp(^-iksQ£lds = 2LeM-^hQ)F{k3Q). 
Substitution of expressions (C-20) and (C-21) into Eq. (C-12) results in 



(C-21) 



p{3) 

fmn 



L^FihQ) 



expf —t(rn — n) j F{k3Q + 7r(?TT- — n)) 



+ exp( i{m — n)]F{k3Q — TT{m — n)) 



TT 



exp( — z(m + n)]F{kzQ + 7r(m + n)) 



vr 



+ exp(-— z(m + n)]F{k3Q - Tr{m + n)) 



In the trigonometric form, this equality reads 

TT{m — n) 



p(3) 
Fmn 



L'F{k3Q){<cos 



cos ■ 



TT{m + n) 



+i< sm 



iT{m — n) 



F{k-iQ + 7r(m - n)) + F{k-iQ - 7r(m - n)) 
Fik^Q + 7r(m + n)) + Fik^Q - 7r(m + n)) 
F{k3Q + 7r(m - n)) - ^(fesQ - TT{m - n)) 



sm ■ 



TT{m + n) 



F(A;3Q + 7r(m + n)) - F{k3Q - 7r(m + n)) 



(C-22) 



Combining Eqs. (C-18) and (C-22), we find that 



p(l) + p{3) 



L 



'F{k3Q){ 



cos ■ 
7r(?7i + n) 



TT{m — n) 



F{k3Q + 7r(m - n)) + ^(/cgQ - 7r(m - n)) 



cos ■ 



F(/c3Q + 7r(m + n)) + ^(fcaQ - 7r(m + n)) 



We now calculate the integrals 



D 



„— , exp( ibi/csQ— ) sin — - — as. 
^ L'^ F 



(C-23) 



(C-24) 



wins 1 



Bearing in mind that sin ^^j^ = g: 



expim 



exp 



we find that 



D^„ 



1 

2i 
L 
2" 



/ exp(i(A;3(5 + vrm) — jds — / exp(z(A;3(5 — vrm) — jds 



1 



k3Q — Tim 
1 



exp(z(A;3(5 ~ vrm)) — 1 
exp(z(A;3(5 + vrm)) — 1 



It follows from this equality and Eq. (C-15) that 



D+ = Lzexpf-A;3(5j expf-— «mJF(A;3(5 - 7rm) - expf — zmjF(/c3Q + 



7rm 
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Similarly, we find that 



Dm = Liexp{--k3Q 



exp( im\F{k'iQ + irm) — exp( —im]F{k^Q — -Km) 



Substitution of these expressions into Eq. (C-12) implies that 



p(2) 



V 



eyip ( —imj F {ksQ + vrm) — exp( im]F{k^Q — irm) 



eyip( --mjF{ksQ + -rm) - exp(-in)F{ksQ 



nn] 



Taking into account that exp( ^m) = cos ^ + isin ^, we arrive at the formula 



p(2) 



L' 



■Km 



cos ( F{k^Q + 7rm) — F{k^Q — nm] 



■Km 



+1 sin — [FiksQ + vrm) + FiksC 



TTm 



vrn 



cos — \F{k^Q + Tin) - F{k3Q - ■ku) 
sin — {F{k:iQ + vrn) + ^(/caQ - im) 



It follows from this equality that 



„(2) 



^ yrnri ~r ISmnjj 



(C-25) 



where 



■Km KTi 

cos cos — 

2 2 

vrm vrn 

+ sin sin — 

2 2 

Km Kn 

sin cos — 

2 2 

TTTiT, vrn 

— cos sm — 

2 2 



As r. 



nm — ' mn 



f^m.n, anQ Sr. 



F{k^Q + 7rm) - F{k^Q - T:m) Fik^Q + vrn) - ^(/csQ - vrn 

FiksQ + 7rm) + ^(/csQ - vrm) ^(fcgQ + to) + FiksQ - ■nn] 
F{k^Q + Km) + F(/c3Q - Km)] \F{k^Q + vrn) - F{k^Q - vrn) 

^(fcaQ + Km) - F{ksQ - -Km) F{ksQ + ttu) + ^(A^sQ - vrn 
Smn) we find from Eq. (C-25) that 



p(2) + (2) 
Irmn ' fnim 



2L^{ 



■Km 



vrn 



cos cos — 

2 2 



-^(A;3(5 + vrm) - ^(/csQ - nm) 



FiksQ + Trn) - Fik^Q - nn] 
nm, . nn 



+ sm sm — 

2 2 



F(A;3Q + vrm) + ^(fcaQ - vrm) 

X FiksQ + vrn) + FiksQ - to)] }. 

Substitution of expressions (C-23) and (C-26) into Eq. (C-11) results in 

vr(?7i — n) 



(C-26) 



Pmn = 2L^[F{k3Q) 



cos ■ 



Fik^Q + vr(m - n)) + F{k^Q - vr(m - n)) 



cos ■ 



vr(?T7. -|- n) 



(F{k^Q + K{m + n)) + Fik^Q - vr(m + n))) 



cos cos 

2 2 



FihQ + -Km) - Fl^ksQ - Km 

Km vrn 



1 KTfl ■Kfl r 

FiksQ + to) - F^ksQ - to)J - sin — sin — [F^k^Q + 
+F{k^Q - vrm)] [F{k^Q + to) + F{k^Q - to)] }. 
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vrm 



(C-27) 



It follows from Eqs. (C-10) and (C-27) that the potential V{r) is given by Eq. (35), where the 
coefficients 

Pmn = l^j U,{k)kUkj p^n{kQx){l - X^)dx (C-28) 

are determined by Eq. (36). 

Setting n = m in Eq. (C-27), we find that 

Pmm = 2L^[F{k^Q)\2F{kzQ)-cos{-Km)(F{k:iQ + 2Tym)+F{k^Q-2T:m) 



■Km 



cos 



Fik^Q + vrm) - Fik^Q - ■nm) 



■Km 



sm 



F{k3Q + TTm) 



+F(A:3Q-7rm)]^}. 
Bearing in mind that 

TTfyi / \ 2 ITTY) / \ 2 

sin^ — [Fik^Q + 7rm) + ^(/cgQ - vrm) j + cos^ — [Fik^Q + 7rm) - ^(/cgQ - Trm) j 
= F'^iksQ + Km) + F'^{k3Q - -Km) - 2 cos{Km)F{k3Q + Km)F{k3Q - -Km), 
we arrive at the formula 

p^m = 2L'^[[2F^{k3Q)-F\k3Q + Km)-F^{k3Q-Km) 

- cos(7rm) F{k3Q) (^(/cgQ + 2Km) + F{k3Q - 27rm)) 

-2F{k3Q + Km)F{k3Q - vrm)] } . (C-29) 



Appendix D 

It follows from Eqs. (43), (C-28) and (C-29) that 

roo rl r °° 1 r 

S = U^{k)k^dk {l-x^)dxlJ2^P^i^Q^ + ^^)+P^i'^Q^-^^ 



m=l 



oo /ixr 

+F{kQx) Y^ ^ 



m=l 



m^ 



F{kQx + 27rm) + F{kQx - 27rm) 



+2 J2 - — ^F{kQx + Km)F{kQx - Km) - 2F^{kQx) ^ 



Bearing in mind that 



m^ 



m=l 



m-^ 



~ J__7r2 



m=l 



m^ 6 



we arrive at Eqs. (44) and (45). 

We begin with calculation of si(z). It follows from Eq. (37) that 



F{z + 2Km) + F{z - 2Km) 



1 . fZ \ 1 . fZ 

sm — |- TTTTi H sm Krn ] . 

z + 2Km V2 J z-2Km \2 



Taking into account that for any integer m. 



(z \ 

— lb vrmj = (— l)'"sin ■ 



(D-1) 
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we obtain 



F{z + 27rm) + F{z - 2-Km) = (-1)™ sin 



Vz + 2: 



1 X (-l)"'z z 

= 2—r — sm — . 

2\z + 27Tm z — 2-Km^ Z'^ — [2-KmY 2 



+ 



This equality together with Eq. (45) imphes that 



.i(z)=2zsiniX:^-,^-^-^^— ^. 



It follows from the identity 



rn^^z"^ — (27r?Ti-)^) z^ ^m? z^ — (27rm) 



+ 



(2vr)2 



that 



sx{z) = -sm- 



OO 1 OO 1 



171=1 



m^ 



■f z"^ — {2Trmy 



— I- 4 > 

6 ^ z2 - (27rm) 

m=l ^ ' 



Using Eqs. (37) and (D-1), we transform this equality as follows: 

siiz) =2ir'^F{z) 
Setting z = 2zi, we obtain 

4E 

m- 

The latter sum is well-known 07j 



(D-2) 



1 



1 



_ = y . 

' z"^ — (27rm)2 , z? — (irm)^ 

1 ^ ' m=l i \ J 



E 

m=l 



1 



z^ — (vrm,)^ 22; 



— ( cot z ) . 

lz\ Z' 



(D-3) 



Substituting expression (D-3) into Eq. (D-2) and returning to the initial notation, we arrive at the 
formula 



1 1 



z 2 
-6 ' zV"2 ~ ~z 
It follows from Eqs. (37) and (45) that the function 82(2;) reads 



Si z =2vr2Fz +- cot--- 
Lh z\ 1 z 



S2iz) = E 



1 psin^(f + ^) sin^(f-^) 



m=l 



m^ 



+ 



(z + 7rm)2 (z — -Km) 



where 



4^^(-) = E 



00 2 I / \2 

z + (irm) 



s^Hz) 



E 



S2 \Z) — $2 [Z) COS Z, 

(-l)™(z2 + (7rm)2) 



''-, m?iz'^ — (vrrTT,)^)^' ^-^ rn^iz"^ — (vrrTT,)^)^ 

m=l ^ \ J / m=l ^ \ J / 

To determine the function Sg yz), we present the first equality in Eq. (D-6) in the form 



(D-4) 



(D-5) 



(D-6) 



s^'\z) 



E 



+ -^E 



' m?[z'^ — (TTm)^)'^ _ (z^ — (vrm)^) 



^2^2- 



(D-7) 
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The first sum in Eq. (D-7) is transformed with the help of the identity 

1 1 



1 ^2 ^2^2 

+ ^ ; TTT + 



oo { -\\m oo f -\\r 

•^2 {Z) - Z 2^ ^2(^2 _(^^\2\2^^ Z^ 



Bearing in mind that 



m=l 

we find from Eq. (D-13) that 



^, m2 \ (2m - 1)2 + \ (2m)2 ' 

m=l ^ ' m=\ ^ ' 



The other sum in Eq. (D-12) is well-known |47j . 



■^ z"^ — (Trmy 2z\smz Z' 
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(D- 



m2(z2 _ (^7]-j72)2)2 ^4 [^^2 ^2 _ (7^772)2 {^£2. _ ^jj-^2'j2 

Combining Eqs. (D-7) and (D-8), we obtain 

S2\z) = ^ E ^ + ^ E -2 7 V2 + 2^^ E T^ 7 vlYl- (D-9) 

m=l m=l ^ ' m.=l ^ \ 1 1 

Differentiation of Eq. (D-3) with respect to z implies that 

E°° 1 _ 1 /cotz 1 2\ m in^ 

_, (z2 - (^m)2)2 ~ 4i2 V^ + ^^z ~ ^)' ^ ^ ^ 

m=l ^ \ / / 

Substitution of expressions (D-1), (D-3) and (D-10) into Eq. (D-9) implies that 

(I) , , vr^ / 1 cot z 1 3 \ ,^ X 

(2) 

According to Eq. (D-6), the function S2 (z) is given by 



-^ m'^iz^ — (7r7?i)2)2 ■^ (2;2 — (7rm)2)2 

771=1 ^ \ y / m,=l ^ \ y / 

Combining this equality with Eq. (D-8), we find that 

s^^\z) = 4 E -^ + ^ E 2 ; ^2 + 2vr' E , 2 / ^2^2 ^ (^-l^) 

z^ ■^ 771^ 2;^ '^ z^ — (Trm)^ '^ (z^ — (vrmi^r 

m=l m=l ^ ' m=l ^ \ / y 

To evaluate the first sum in Eq. (D-12), we present Eq. (D-1) in the form 

200 -I 00-1 00 1 -lOO-i 

Y ~ ^, (2m - 1)2 ^ ^^ (2m)2 ~ ^_, (2m - 1)2 "^ 4 ^, ^' 

m=l *■ ' m=l ^ ' m=l ^ ' m=\ 

Combination of this equality with Eq. (D-1) results in 

^2 00 -. 2 



^, (2m - 1)2 ~ y ^_, (2^ ~ 24' ^ ^-^^^ 

m=l ^ ' m=l ^ ' 



V ^^ = --. (D-14) 

^ m2 12 ^ ' 

rra=l 



■■^-^^ ^2 — {TimX^ 2zVsinz 2/ 



Differentiation of Eq. (D-15) with respect to z results in 



E 



.\y 



1/1 cosz 2 

+ 



^i 



' (z^ — (7rm)2)2 4^2 V 2; sin z sin^ z z^ ) 

m=l ^ \ / y 

Substitution of expressions (D-14) to (D-16) into Eq. (D-12) implies that 



J2), 



cosz 3 



_ TT^ / 1 
■^2 'l^j - 7Tl2(~c ' „„;„ „ ' I7:X, ^27- 



2^2 V 6 z sin z sm" z z^ 



Combining this equality with Eqs. (D-5) and (D-11), we arrive at the formula 

vr^ / 8 + cos z 1 — cos z ^ 

S2{Z) 



2z^\ 6 



It follows from Eqs. (37) and (45) that 



- (-l)- sin(f + ^)sin(f-^ 

^3{Z) — 2_^ 2 



m=l 



m^ 



z"^ — (irm)'^ 



Bearing in mind that 



z Tim z irm s 1 / . n \ 1 

sm( — I sm = - cos 7rm — cosz = - 

^2 2 2 2 2 V ^ ^ / 2 



-D^-cosz 



we find that 



where 



.(1) 



(2), 



,(1) 



S^ '[Z] 



E 



s^iz) = 7; 4 (z) -4 (^) COS 2; 



'-j^ m'^{z'^ — (Trm)^) 



4^^(^) = E 



' m?{z'^ — {-nm)'^) 



Taking into account that 



1 



1 



+ 



vr 



im?{z'^ — (Trm)^) z^ ^m? z"^ — {irm)^ 
and employing Eqs. (D-1) and (D-3), we obtain 

^2.1 1, 1 



vr 



4"(^) = ^ 



- + - cot z 

3 z V z 



Similarly, we find that 



s^'\z 



) = ;2iE^7^ + ^'E 



z^ ^' — \ m^ 

m=l 



■f z^ — (7rm)2- 

Substitution of expressions (D-14) and (D-15) into Eq. (D-20) implies that 

-2^1 1/1 1^ 






+ 



6 z^sinz Z' 
Combining this equality with Eqs. (D-18) and (D-19), we find that 



S3 



«-£.( 



TT /2 + COS z 1 — cos z 



4^2 V 6 
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(D-16) 



(D-17) 



(D-18) 



(D-19) 



(D-20) 



(D-21) 



Substitution of Eqs. (37), (D-4), (D-17) and (D-21) into Eq. (45) results in Eq. (46). 
It follows from Eqs. (9) and (47) that 



S 



vr 



1-x' 



62Q2 Jq X 



dx / a{kbQx)k dk. 







Introducing the new variable z = kbQx, we obtain 



S 



b^Q5 Jo X 



-dx 







a{z)z dz. 



We now set y = k^:bQx and change the order of integration, 



S 



vr 



2^4 i-k. 



bQ Jo 



y \^\dy 



KbQ^ ' r ^0 



bQ Jo 



a{z)z dz 



k*bQ 



1 



a{z)z dz 
V \^\dy 



(D-22) 



Substituting the expression 



k,bQ 



1 



y \^\dy 1/1 



k,bQ^ h^ 4Vz2 (fc,6Q)2 



into Eq. (D-22), we arrive at Eq. (48). 

To transform expression (49) for the function Aq{x), we, first, find from Eq. (46) that 



a{z) = l + 2z^ — 



d /I — cos z 



dz V z^ 
Substitution of Eq. (D-23) into Eq. (49) and integration by parts imply that 



(D-23) 



1 



Ao{x) = -T 



\x^ _ ,^)L^^dz - A I [x'-z^)dz 



Using Eq. (18) and calculating the second integral in this equality, we obtain 



Anix) = -^ 



A{x) 



1 



(D-24) 



Using Eq. (D-24) and expanding the function A{x) into the Taylor series in the vicinity of the 
point a; = 0, we find that 

limAo(x) = 4^(0). (D-25) 

x^o dx^ 

Differentiation of Eq. (A-8) with respect to x implies that 



d^A A /"^ 

dx'^ x^ 7o z^ 



z — sin z , — 8x + 10 sin x — 2x cos x 

dz H = . 



Applying L'Hospital's rule, we obtain 



d^A 
dx"^ 



(0) 



1 

90" 



(D-26) 



The first equality in Eq. (50) follows from Eqs. (D-25) and (D-26). The other equality is a 
consequence of Eqs. (19) and (D-24). 



29 



Appendix E 

Substitution of expression (D-24) into Eq. (53) implies that 



C2 = 86^ 



X. ( 3x^\ 1 rhx 

^(^)^"K-— J^^-eVT 



Our aim now is to prove that the first term in the square brackets 






is small compared with unity. It follows from Eq. (A-8) that for any y, 



Beaing in mind that 

we find from Eq. (E-3) that 



2 /""^ 2 /""^ 

^(y) / a{z)dz = — / a{z)dz. 

y Jo y Jy 



1 / sinzx 2 
0<a(z) = -(l <- 



2 /""^ 
A{y) - - / aiz)dz 

y Jo 



< 



A r'^ dz 4 



y Jy z^ y^ ' 



Setting y = x/6 in Eq. (E-2), we obtain 



roo , QX2„.2 
l = bj ^(y)exp( -^]dy = 5(/i+/2 + /3) 



where 



h 
h 



■,„)exp(-^).. 



U'y 



2 

352y2s dy 



2„,2, 



c^y, 



/•OO /"OO 

Is = 2 a{z)dz I exp 
It follows from Eqs. (19) and (E-6) that 



2 / y 



\h\ < I A{y)dy < i 
JO 3 



Equations (E-4) and (E-6) imply that 



. f°° / "^5 y \dy . f^ ~„ 
\h\<^j^ «^p(-^)i^4/ -1 = 4. 



2 Jy 



dy 
1 y^ 



According to Eqs. (B-14) and (E-6), 



'^=2 A '"P 



3S^y^\dy it f°° / 3x^\dx 
exp 



2 ^ y 2Js 



2 / X ' 



(E-1) 



(E-2) 



(E-3) 



(E-4) 



(E-5) 



(E-6) 



(E-7) 



(E- 
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where we return to the variable x. For 5 < 1, this equahty reads 



'3 = |(4" + /f). 



where 



/. 



(1) 



^ / Sx^ \ dx 
exp 



r(2) 



It follows from these equalities that 









dx 



In (5, 



l/fl<rexp(-^)d.< 



exp ^ I dx 



Is x " Ji V 2 / Jo V 2 

Combining Eq. (E-5) with Eqs. (E-7) to (E-10), we arrive at the formula 



3+21V6-'"* 



which completes the proof. 

It follows from Eqs. (57) and (D-24) that 



B2 



^(T)exp( — — jx^dx 



1 

3 JO 



°^ / 3x 
exp 



{-^)x'dx^ 



(E-9) 



(E-10) 



(E-11) 



(E-12) 



According to Eq. (25), the first integral in Eq. (E-12) coincides with Ci. The other integral is 
given by Eq. (B-1). Substitution of these expressions into Eq. (E-12) results in 



Bo 



Ci 



1 2tt 



18 V 3 



Combining this equality with Eq. (26) and disregarding terms beyond the second order of smallness, 
we arrive at Eq. (58). 



Appendix F 

To find Ai, we formally expand the denominator in the first equality in Eq. (64) into the Taylor 
series in ^i = — ^u, 



Ai[y) 



Ai (y) \l + fiiAi (y) + J2 (/"i^i iv) 



l-fiiAi{y) - ^^2 

and evaluate the integrals of appropriate terms in this sum. The integral 

3x2 



(F-1) 



^i^=fM^)ew{-^ydx 



was calculated previously. It follows from Eqs. (53), (54) and (62) that the leading term in the 
expression for A^ is given by 



A 



(1) _ 2 /2^^2^ 



27 V 3 



(F-2) 



To find the coefficient 



A? = f Al{^)e.p{-^-f)x'dx, 
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we substitute expressions (62) and (D-24) into this equality and obtain 



A 



(2) 



9 roo 

^ ir2 



il* io ^^1- 



X, r . ,x 



■^fj' 



expl ^ \ax. 



(F-3) 



The function 



a\{x) = Aq{x) A{x) 



is even and continuous, it has the zero hmit at x ^ and vanishes being proportional to x^^ at 
X — > oo. This means that the Fourier transform ai(s) of this function exists. Presenting Eq. (F-3) 
in the form 

exp( ;^ "X / Qi(s)expl 



A 



(2) 



1627r ]- 



ISX 



ds 



and evaluating the integral with the use of Eq. (B-7), we arrive at 

X3 /.oo , -,2, 



A 



(2) 



81V27r J- 



ai ((5z v3) exp ( \dz. 

-oo ^ 2 ' 



The leading term in this expression is obtained when the function ai{s) is expanded into the Taylor 
series in the vicinity of s = and only the first term in the series is taken into account, 



(o\ 5'^ 25"^ r°^ r 1 



dx. 



(F-4) 



The integrals 



A 



(m) 



AT{^)e^p{-^-^]x'dx 



with 771 > 2 are estimated by using a standard approach. Setting 

Q„(x) = A^(X), 

we re-write this equality in the form 

3x2. 



aM r (X ( 3x-x 2 , 

A\ ' = Q„(-)exp(^-— jx dx. 



(F-5) 



(F-6) 



It follows from Eqs. (62) and (D-24) that A'i'{x) vanishes at x ^ oo being proportional to x^^™. 
This means that for any m > 2, the second derivative with respect to s exists of the Fourier 
transform 

o.m{s) = / Qm(x) exp(isx)(ix (F-7) 



of the function amix). Bearing in mind that the function amix) is even, substituting the expression 







2tt 



am(s)exp( —)ds 





into Eq. (F-6) and changing the order of integration, we find that 



Al' 



(m) 



Att 



am{s)ds 



exp 



3x isx\ 



X dx. 
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Using Eq. (B-19) for the internal integral and setting z = s/((5v3), we arrive at 



A 



(™) = ^ y_ ""^(^^^) e-p(-y) (1 - ^')rf^- 



(F- 



Expanding the even function OLmi^s) into the Taylor series in the vicinity of the point s = and 
neglecting terms beyond the second order of smallness, we find that 

We combine this equality with Eq. (F-8), use Eq. (B-21), and obtain 



(m) 



(53 d^ 



Ur. 



It follows from Eq. (F-7) that 
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This equality together with Eqs. (F-5) and (F-9) yields 
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Substitution of expressions (F-1), (F-2), (F-4) and (F-10) into Eq. (64) results in 

2 /97r 9 C'^ / 1 \ /) C'^ °^ 

^ ■m=3 

Taking into account that 

and returning to the initial notation, we find that 



A, = 6^ 



2 i2tt 4 f^ / 1 \ f^ 
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Combining the integral terms and utilizing Eqs. (62) and (D-24), we obtain 
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Ai = 6^ 
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(F-12) 



As all integrals in Eq. (F-12) converge, it can be shown that the above transformations are correct 
for an arbitrary (not necessary small) ^. 

Equation (F-12) implies that the leading term in the expression for Ai is (at least) of order 
of (5^. As the quantity A2 is included into Eq. (63) with the pre- factor S'^, in order to disregard 
this term it suffices to show that A2 is small compared with unity. Taking into account that the 
function Ai{x) is non-negative, we write 



A{x) 



l + fj.Ai{x) 
33 



< A{x) 



(F-13) 



This equality together with Eqs. (25) and (64) irtiphes that 

3x2 ■ 



< A2 < I" A{^)e^p[-^-)x^dx = Ci. 



It follows from this inequality and Eq. (26) that the contribution of A2 into Eq. (63) is negligible 
compared with Ai. This conclusion together with Eqs. (63) and (F-12) implies that 
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Equation (65) follows from this equality and Eq. (62). 

To find the coefficient Ei, we evaluate the integrals of appropriate terms in expansion (F-1). 
The integral 

r^^^ = f A,Qe^p{-^-^)x'dx 

was calculated previously. It follows from Eqs. (57), (58) and (62) that the leading term in the 
expression for T\ reads 
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Using Eqs. (62) and (D-24), we transform the coefficient 
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as follows: 
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The first term in the square brackets is of order of 6, because the Fourier transform of the function 
A'^{x) exists. The second term is estimated in Eq. (E-11), where it is shown that it is of order 
of (5|ln(5|. The last term is calculated explicitly, and it is of order of unity. Neglecting small 



(2) 
contributions into the coefficient E^ , we obtain 
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For any integer m >3, the coefficient 



r\ 



(m) 



'AT{^)exp[-^]x^dx 



reads 



j.(™) ^ ^4 



/ /3m(^)exp( — ^)dx, 



where Pm{x) = ^™(x)x^. As A^^^x) is an even continuous function, and it decreases being propor- 
tional to x~2"^ as X ^ 00, the Fourier transform j3m{s) exists of the function /3m (x). Evaluating 
the integral by analogy with Eq. (B-3), we find the leading term in the expression for r^™\ 
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It follows from Eqs. (F-1) and (F-14) to (F-16) that 
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Combining this equality with Eq. (F-11) and returning to the initial notation, we arrive at the 
formula 
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To assess the coefficient r2, we use Eqs. (25), (26) and (F-13) and obtain (up to terms of higher 
order of smallness) 

o<r.<l Ai-){-—)x^dx = -, 

which implies that the contribution of r2 into the mean square end-to-end distance is negligible 
compared with that of Fi [the integral F2 is multiplied by 6'^ in Eq. (66)]. Substitution of Eq. 
(F-17) into Eq. (66) results in 
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Combining this equality with Eq. (62), we conclude that the second term in the square brackets 
is negligible compared with the first. Neglecting this contribution and returning to the initial 
notation, we arrive at Eq. (68). 

To assess the integral terms in Eqs. (65) and (69), we introduce the function 

Al{x)x'^ 
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l + iJ.Ai{x) 



dx. 
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If the function Ai(x)x'^ had been integrable, we could use the inequality 



fJ-Aiix) 



< 1 



to find that D(fi) < /q°° Ai(x)x'^dx < 00 for any positive fi. As this is not the case, more sophisti- 
cated estimates are needed. It follows from Eqs. (62) and (D-24) that 
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Substitution of this expression into Eq. (F-19) yields 
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As the function A{x) monotonically decreases with x and vanishes at x ^ 00, there exists an xq 
such that 

I < Aix) < ^ (0 < X < xo). 



< Aix) < - (xo < X < 00) 
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We present the integral in Eq. (F-20) as the sum of two integrals: from zero to xq, and from xq to 
infinity, and evaluate appropriate integrals separately, 
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Taking into account that in the interval x £ [0,xo], 

x' + 1^(1 - Aix)) > !/.(! - Aix)) , I - Aix) < i, 

the function Z?i(/u) is estimated as foUows: 

9 f^o / 1 \ 3 

D,it,)<-J^ ^--Aix))dx<-xo. (F-23) 

Using the estimates in the interval x £ (xq, oo), 

x^ + -n( A(x)] > x^ H a, A(x) < -, 

the function D2{^) is evaluated as 

^'^^^-9X0 x2 + ^^-9Jo x^ + §^fi^9^2^Jo ^^Tl = 2V6- ^^"^^^ 

It follows from Eqs. (F-22) to (F-24) that there are positive constants D^^' and D^^' such that for 
any ^, 



(5V27r 

where notation (62) is employed. This formula implies that the integral term in Eqs. (65) and (69) 
obeys the inequality 

e /•- _A^ixW^^ ^ 27r<5^ H -^^^dx < 2^Z)(o)<5 + 15(1)^2^, (F-25) 

Jo l + /uAi(x) Jo l + iJ.Ai{x) 

and it can be disregarded compared with unity. 
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Figure 1: Graphs of the functions A{x) and Aq{x). 
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